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It is investigated under which conditions an adiabatic adaption of the dynamic and spectral 
information of vector mesons to the changing medium in heavy ion collisions, as assumed in schematic 
model calculations and microscopic transport simulations, is a valid assumption. Therefore time 
dependent medium modifications of low mass vector mesons are studied within a nonequilibrium 
quantum field theoretical description. Timescales for the adaption of the spectral properties are 
given and non-equilibrium dilepton yields are calculated, leading to the result that memory effects 
are not negligible for most scenarios. 



I. INTRODUCTION AND MOTIVATION 

High energy heavy ion reactions allow for studying strongly interacting matter under extreme conditions, i.e., 
high densities and temperatures. Photons and dileptons do not undergo strong interactions and thus may carry 
undistorted information on the early hot and dense phases of the fireball, because the production rates increase 
rapidly with temperature. Dilepton spectra are expected to play a central role in inferring the restoration of the 
spontaneous breaking of chiral symmetry from heavy ion reactions. In the low mass region they couple directly to the 
light vector mesons and reflect their mass distribution. They are thus considered the prime observable in studying 
mass (de-)generation related to the restoration of spontaneous chiral symmetry breaking. 

In this work we investigate dilepton production from the hot hadronic medium and we will concentrate on the low 
mass region, particularly the medium modifications of the p-meson. The medium and hence also possibly the properties 
of the regarded mesons undergo substantial changes over time. Such scenarios have been described within transport 
calculations using some quantum mechanically inspired off-shell propagation 0,0- It emerges the important question, 
whether a quasi instantaneous adaption of the dynamic and spectral information to the changing medium, as assumed 
in more schematic fireball model calculations |3| and microscopic transport simulations, is a suitable assumption or 
whether the vector meson's spectral information reacts to changes with a certain " quantum mechanical" retardation. 
The transport description of off-shell excitations is an open field of research and necessary in order to understand the 
transport dynamics of resonances. We employ a nonequilibrium quantum field theoretical description based on the 
formalism established by Schwinger and Keldysh 0, Q . We give a formula for the dynamic dilepton production rate 
and simulate modifications of the light vector mesons due to the dynamically changing medium in heavy ion collisions 
by parameterizing a certain time dependence of the p-meson self energies. We are able to analyze the mesons' dynamic 
spectral properties as well as the resulting dilepton rate and the yield from an evolving fireball and compare to the 
quantities computed assuming adiabaticity. 



II. THE NONEQUILIBRIUM PRODUCTION RATE 



We utilize the Schwinger-Keldysh formalism in order to derive the dynamic non-equilibrium rate of produced 
electron-positron pairs, coming from the decay of light vector mesons via virtual photons in a spatially homogeneous 
system (details in 0). 

Projecting on the particle number in the electron propagator G < and using the equations of motion for G < , the 
so called Kadanoff-Baym equations, we find the production rate of electrons for a homogeneous, yet time dependent 
system, to read 
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with the electron self energy S and po = E p . The free electron propagator can be used because due the electrons' 
long mean free path the electrons are not expected to interact with the medium after being produced. 

The medium effects enter via the dressing of the virtual photon propagator in the electron self energy (see Fig. 
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FIG. 1: £(1,2) in coordinate space 



D-y is the virtual photon propagator and II is its self energy. We have 
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with k the momentum of the virtual photon. On inserting this self energy and defining p + — k — p as the four- 
momentum of the positron and p~ = p as that of the electron, equation Q becomes 
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Applying the equilibrium properties of D < , it can be shown that equation © is the generalization of the well known 
thermal production rate for lepton pairs in the stationary case : 
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In the following, we will consider the mode k = exclusively, i.e., the virtual photon resting with respect to the 
medium. After projecting on the virtual photon momentum and taking the electron mass to zero, we get 
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The dynamic information is inherent in the memory integral on the right that runs over all virtual photon occupation 
numbers D< from the initial time to the present. This way the full nonequilibrium electron production rate at 
the present time r is determined. We introduce the dynamic medium dependence by dressing the virtual photon 
propagator with the medium dependent p- or w-meson. This dressing enters with the self energy II < via the relation 
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that follows from the Kadanoff-Baym equations. Using projectors on the different polarizations and the equality of 
each polarization (k = 0), we find 
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For dilepton production II rct oc e 2 and E is the invariant mass of the virtual photon. For the cases we are interested in, 
|II ret | <C E and we can approximate = Dg et (8)11^ <S> Dq v . Diverging contributions at early times (low frequencies), 
due to the undamped photon propagators, cause numerical problems. We introduce an additional cutoff A for these 
propagators: Dq*(t — t\) = (r — t\) — » (r — ti)e _A ( T_tl ) and analogously for {t 2 — t). The exponential factors 
lead to a reduction of the rate, which we will overcome by renormalizing the final result. This does not affect the 
timescales we are interested in, and comparison of the dynamically computed rate for a stationary situation (constant 
self energy) with the analytic, thermal rate shows perfect agreement. 
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Vector meson dominance (VMD) allows for the calculation of the photon polarization tensor ILj,, using the identity 
between the electromagnetic current and the canonical interpolating fields of the vector mesons, which leads to 
= p-TOp-D^ a /3- Once more, we apply the generalized fluctuation dissipation relation © to calculate 

D< T = D* p %®Z< T ®D™ T , (8) 

with the p-meson self energy £<. The transverse parts of the retarded and advanced propagators Z?™^(k, t\, t$) = 
_Dp V T (k, ts, tx) of the vector meson in a spatially homogeneous and isotropic medium follow the equation of motion 

(-3* - m 2 - k 2 ) D p %(k, t u t 3 ) - f 1 diZf T (k, h, t)D*$(k, t, h) = 5{t x - h). (9) 

Jt 3 

The dynamic medium evolution is now introduced by hand via a specified time dependent retarded meson self energy 
vjrct(i-,w) ||| From that the self energy S < , needed for solving equation (JSJ), follows by introducing a background tem- 
perature of the fireball. The fireball is assumed to generate the time dependent self energy S rot and, assuming a nearly 
quasi thermalized system, the p-meson current-current correlator E < is given via S < (r, to) — 2inB{T{T))ImY7 ct {T,uj) 
, which follows from the KMS relation ||, being valid for thermal systems. The latter is a rather strong assumption, 
but necessary in order to proceed. 



III. THE MEDIUM IN NONEQUILIBRIUM 

The medium effects are introduced via a specific evolving self energy of the vector meson. A simple self energy 

£ rct (w,r) = -iojT(r), (10) 

with a k- and w-independent width T, leads to a Breit-Wigner distribution for the spectral function. The time 
dependence is being accounted for by introduction of the parameter r. For the k = mode, the full self energy for 
coupling to J p = i -resonances is given by 
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with E = m 2 R + k 2 — tojv and rriR and mjy the masses of the resonance and the nucleon respectively. Tr is the 
width of the resonance and gi the isospin factor . For our purpose it suffices to retain the part of the self energy 
that creates the pole structure. We will neglect the uj 2 in the numerator because it causes straightforward dispersion 

relations to become invalid (subtracted dispersion relations are needed in this case). We will shorten ^™ P ^J j;/2 w 2 , 
with uj 2 = lGeV 2 by C, a dimensionless factor. 




FIG. 2: Linear switching off of in-medium effects over a certain time At 
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Investigating which contributions in equation come from which times in the past, shows that there are contri- 
butions from early times as well as alternating positive and negative contributions. An interpretation of this becomes 
difficult and it follows that only the time integrated yield is a physical quantity. 

In order to quantify the times that the mesons' spectral properties need to adjust to the evolving medium we 
change the self energy linearly in time (see Fig. |5J). As a possible characteristic timescale we consider the difference 
of the final spectral function to the dynamically calculated one at the time where the medium effects arc fully turned 
off, described by the difference in the moment J Q A{w) 2 u 2 du> of the two spectral functions or the difference in the 
peak position and height. All methods lead to similar results 6]. We find an exponentially decreasing difference with 
increasing duration of the change At (see Fig. |2J), from that we extract a time constant f. For the the p-meson we 
find a typical timescale of about 3 fm/c. That means that the spectral properties follow the changes in the medium 
nearly instantaneously only if the evolution is very slow as compared to 3 fm/c. This means that the vector mesons 
possess a certain memory of the past, and even if they decay outside the medium, they still carry information on 
the medium in that they were produced. This becomes important especially for the w-meson, having a width of only 
8.49 MeV (f w 60 fm/c). It turns out that f is proportional to C/T2 with c lying between 2 and 3.5, depending on 
m and Ti. This is significantly longer than the naively expected timescale of 1/r. The time needed by the dilepton 
rate to follow changes is approximately equal to that of the spectral function. The quantum mechanical nature of the 
regarded systems leads to oscillations and negative values in the changing spectral functions, occupation numbers and 
production rates as well as interferences that one does not get in semi-classical, adiabatic calculations. The "rate" 
calculated here possesses the full quantum mechanical information incorporated and contains " memory" interferences 
that might cause cancellations - hence the rate has to be able to become negative while the time integrated yield 
always stays positive as the only observable physical quantity. An intriguing example for the occurring oscillations is 
shown in Fig. [21 
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FIG. 3: Production rate for the change of the mass from m = 400 MeV to 770 MeV (constant T=150 MeV and constant 
T — 160 MeV) directly after the self energy has reached its final form (after 7.8 fm/c) and 2 (and 4) fm/c later. Oscillations 
and negative values appear in the intermediate rates (left). The corresponding yield stays positive (right). 



To calculate the yield, we model the fireball evolution and fold it with the calculated time dependent rates, similar 
to 3J. For the effective volume we choose a longitudinal Bjorken expansion together with an accelerating radial flow 

V eff (t > t ) = irct{r + v (t - t ) + Q.5a (t - t ) 2 ) 2 , (12) 

with ro = 6.5 fm, vq = 0.15 c and ao = 0.05 c 2 /fm (see also ^3)- From <I12|) and the constraint of conserved 
entropy (given by a constant entropy per baryon S/A — 30 for SPS energies [Tl|), temperature T(t) and chemical 
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potentials follow as functions of time. We start the calculation at the freezeout temperature of 175 MeV, whereas 
the final temperature, reached after a lifetime of about 7.8 fm/c is 120 MeV (thermal freezeout). At this point, we 
turn off further dilepton production by a rapid decrease of the temperature towards zero. With the time dependent 
temperatures we can integrate the rate and immediately find the yield per unit four momentum. The results for 
different scenarios are shown in Fig. We compare to Markov calculations that assume instantaneous adaption of 
the spectral function (and the rate) to the self energy, as employed in Q (dashed lines in FigQJ. The yield resulting 
from a mass shift following from assumed Brown-Rho scaling (l2l | shows an enhancement of about a factor of 2 in the 
dynamic calculation within the mass range where the CERES experiment |13l Il4j has measured a strongly enhanced 
dilepton yield compared to calculations with vacuum spectral functions. Also the coupling to the N(1520) resonance 
with no broadening shows an enhanced production around the resonance peak and the p vacuum peak, but also a 
reduced yield in the region between the peaks. On the other hand, the difference to Markov calculations becomes 
smaller for large in-medium widths. 




FIG. 4: Comparison of the dynamically computed dilepton yield from the fireball (solid line) to the one calculated assuming 
adiabaticty (as done in principle in 3]) (dashed line), to show the differences caused by memory effects. Mass shift to 400 MeV 
in-medium mass (left) and coupling to the N(1520) resonance (right). No broadening of either the resonance or the p-meson is 
included here. Also shown is the yield coming from a constant p meson's vacuum spectral function (dotted line) for comparison. 



IV. SUMMARY AND CONCLUSIONS 



In the present work we introduced a method to calculate dilepton production rates within a non-equilibrium field 
theory formalism, based on the real time approach of Schwinger and Keldysh. We investigated possible medium 
modifications of the p meson in a fireball created in a heavy ion collision. We considered mass shifts, broadening and 
coupling to resonances. Special attention was put to possible retardation effects concerning the off-shell evolution. 

The timescale on that the spectral function adjusts to changes in the self energy was found to be proportional to 
the inverse vacuum width of the meson T 2 like c/r 2 , with c approximately 3. Further dependence on the in-medium 
width as well as on the size of the medium modification is present. We find typical retardations for the p of 3 fm/c 
and about 60 fm/c for the w, a very slow adjustment. 

The full quantum field theoretical treatment leads to oscillations in all mentioned quantities when changes in the 
self energy are performed. This oscillatory behavior reveals the quantum mechanical character of the many particle 
system, present in the investigated heavy ion reaction. The oscillations potentially cancel when the rate is integrated 
over time such that the measurable dilepton yield is always positive. 
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Comparison of dynamically calculated yields with those calculated assuming adiabaticity reveals differences. 
About a factor of 2 difference was found within the invariant mass range of 250 to 500 MeV for mass shifts predicted 
using Brown- Rho scaling for the p-meson in a fireball at SpS energies (158 AGeV). This is the range where CERES 
measures an increased dilepton yield as compared to calculations assuming the p's vacuum shape. Similar results 
were found for the coupling of the /9-meson to resonance-hole pairs. Our findings show that exact treatment of 
medium modifications in principle requires the consideration of memory effects. Further investigation of the w-meson 
is being done - due to its small width it causes numerical complications. It has to be seen how our results can be 
consistently (and probably in an approximative manner) incorporated in semi-classical nonequilibrium transport 
codes like UrQMD [H, BUU Q or HSD pjj. 
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